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Abstract. In analogy with ocean waves running up towards the beach, shoaling of 
prechirped optical pulses may occur in the normal group-velocity dispersion regime of 
optical fibers. We present exact Riemann wave solutions of the optical shallow water 
equations and show that they agree remarkably well with the numerical solutions of 
the nonlinear Schrodinger equation, at least up to the point where a vertical pulse 
front develops. We also reveal that extreme wave events or optical tsunamis may be 
generated in dispersion tapered fibers in the presence of higher-order dispersion. 
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1. Introduction 

Extreme or rogue waves have received a great deal of attention recently for their 
emergence in a variety of applications, ranging from fluid dynamics and oceanography 
to plasma physics, Bose-Einstein condensation (BEC) and nonlinear optics [I]-[2]- The 
most popular manifestation of a rogue wave has so far been described as the sudden 
build-up and subsequent rapid disapperance in the open sea of an isolated giant wave, 
whose height and steepness are much larger than the corresponding average values of 
other waves in the ocean. A different manifestation of rogue waves with the potential 
for large scale damages also occurs as a result of wave propagation in shallow waters, 
a framework which describes the run-up of tsunamis towards the coast [ L 3J - Moreover, 
it is also known that the crossing of waters propagating in different directions or with 
opposite velocities may lead to the formation of high-elevation and steep humps of 
water or sneaker waves, which may cause severe disruptions along the coastine and river 
flooding jl]. The same phenomena also occur in the run-down of avalanches falling from 
a mountain channel or a glacier [5] . 

In deep waters, the dynamics of rogue wave formation may be described by the well- 
known one-dimensional nonlinear Schrodinger equation (NLSE). As such, the generation 
of rogue waves has been closely associated so far with the presence of CW breaking or 
modulation instability (MI) |6j, which occurs in the so-called self-focusing or anomalous 
group-velocity dispersion (GVD) regime. Thanks to the integrability of the NLSE, 
the nonlinear development of the MI originates families of exact solutions such as the 
Akhmediev breathers [7\. Although generally periodic both in the evolution variable 
(e.g., distance) and the transverse dimension (e.g., time), such nonlinear wave families 
may also include strictly spatio-temporally localized waves: consider for example the 
Peregrine soliton [8]. 

Given the widespread applicability of the NLSE to fields of physics other than 
oceanography, it has been possible to predict and experimentally observe rogue wave 
phenomena in different contexts such as nonlinear optical fibers. For example, the 
temporal statistics of optical supercontinuum generation has revealed the emergence of 
extreme solitary wave emissions [S]. In addition, optical fibers have provided the test- 
bed for the the first clear experimental observation of Peregrine solitons by inducing the 
MI of a pump laser by means of an additional seed laser [10J. 

On the other hand, the appearance of rogue waves in shallow waters is a field 
which is still the emerging stage [3]. As such, it does not appear to have received 
much attention in other physical contexts besides hydrodynamics. A recent study has 
however pointed out that extreme waves may be also be generated in optical fibers in 
the normal GVD regime of pulse propagation, where MI is absent [TTJ . In this famework 
the model linking hydrodynamics with nonlinear optics is provided by the semiclassical 
approximation to the NLSE [12], which is known as the nonlinear shallow water equation 
(NSWE) [T3] or the Saint- Venant equation [2]. Although in the normal GVD regime a 
CW is modulationally stable, extreme waves may still be generated in the presence of a 
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suitable temporal pre-chirp or phase modulation |llj.|15]-|16j. which is analogous to a 
nonuniform velocity distribution of the propagating water waves, eventually leading to 
tsunamis. As well known, the conditions of shallow water (as opposed to deep water) 
propagation applies in oceanography whenever the wavelength of the wave is much 
longer than the depth of the water, such as it happens for tsunamis even at large 
distances from the coast, given their typical wavelengths of the order of thousands of 
km. On the other hand in nonlinear optics, the applicability of the NSWE requires the 
characteristic dispersion distance to be much longer than the nonlinear distance (small 
dispersion limit) [12]. 

In this work we show that, in analogy with the commonly experienced case of 
ocean waves as they run-up to the beach, the shoaling of properly pre-chirped optical 
pulses may also occur in the normal dispersion regime of optical fibers. Therefore we 
shall consider the propagation of special, input temporally pre-chirped optical pulses 
with different power profiles. These pulses represent nonlinear invariant solutions of the 
NSWE (Riemann waves). For such type of chirped pulses, we obtain exact solutions 
of the optical NSWE, and demonstrate their good agreement with numerical solutions 
of the NLSE, at least up to the point where a vertical pulse edge or front develops 
in the power profile. In addition, we also present simulations which reveal that third- 
order dispersion (TOD) leads to the occurrence of extreme waves or optical tsunamis 
whenever a dispersion tapered fiber is used, in analogy with the dramatic run-up and 
wave height amplification of a tsunami as the coast is approached and the water depth 
is progressively reduced. 

In section [2] we present the NSWE and its formulation in Riemann invariant form, 
which leads to the definition of the nonlinear Riemann waves. Exact solutions for these 
temporally chirped optical Riemann pulses with different input power temporal profiles, 
such as the parabolic, Gaussian or hyperbolic secant profile, are presented in section 
[3| and critically compared with the numerical solutions of the NLSE. The analysis of 
section [3] reveals the emergence of wave-breaking free optical wave shoaling, that is the 
formation of a vertical edge for the propagating pulse, without any appearance of high 
frequency temporal oscillations or shock phenomena. Finally, as discussed in section 
[4j we show that whenever the balance of GVD and TOD dynamically evolves along 
the fiber, such as it occurs in dispersion varying or tapered fibers, significant temporal 
compression and peak power amplification may result in the form of a high-intensity 
spatio-temporally localized flash of light propagating entirely in the normal GVD regime. 

2. Basic equations 

As it is well known, the propagation of short light pulses in optical fibers may be 
described in terms of the NLSE 



.dQ_f32(PQ_ d 3 Q 
1 dz 2 dt 2 \ dt 3 



+ l\Q\ 2 Q = ctftQ. 



(1) 
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Here z and t denote the distance and retarded time (in the frame travelling at the 
carrier frequency group- velocity) coordinates; (3 2 , (3^ and 7 represent GVD, TOD and 
the nonlinear coefficient, respectively; Q is the field envelope. Moreover, atf denotes the 
rate of the carrier frequency shift with distance. A continuous frequency shift may be 
introduced in a fiber ring laser by means of an acousto-optic filter [12]. Otherwise, the 
frequency shifting term may be introduced to represent, in combination with the TOD 
term, propagation in a dispersion tapered optical fiber, where the local dispersion varies 
along the propagation distance [18]- [19]. In dimensionless units, and in the normal GVD 
regime (i.e., (3 2 > 0), Eq.Q reads as 



dq (3 2 dq .{3 d 3 q 



i — ^ — + \q\ Q — oTq. 



(2) 



dZ 2 dT 2 6 dT 3 
Here T = t/t , Z = z^P = z/L NL , (3 2 = /3 2 /(*o7-fo) = L NL /L D , where L NL and L D 
are the nonlinear and dispersion lengths, respectively, (3 = (3%/ (t 3 ) 'jPo), q = Q/yfPo, and 
a = tt/t /(7-Po); ^0 an d -Po are arbitrary time and power units. Eq.([2| can be expressed 
in terms of the real variables p and u which denote the field dimensionless power and 
instantaneous frequency (or chirp) 



q( T , z ) = \M T > Z)exp 



u(T', Z)dT 



(3) 



By ignoring higher order time derivatives in the resulting equations (which is 
justified for relatively small values of (3) as well as TOD, one obtains from the NLSE 
the hydrodynamic NSWEs 



(4) 



dp d (pu) 



dZ' 

du 
dZ> 



+ 



+ u 



dT 

du 
df 



= 

dp 
df 



a. 



where Z' = (3Z and a = a/(3. In hydrodynamics, Eqs.Q describe the motion of 
a surface wave in shallow water, i.e., a wave whose wavelength is much larger than 
the water depth. In this context, p and u represent the water depth and its velocity, 
respectively. Moreover a = dh/dT, where h{T) is the unperturbed water depth along 
the channel or beach axis [3]. Therefore a constant a represents tsunami run-up towards 
a beach with uniform slope. 

Solutions of Eqs. Q may be found by making use of the Riemann invariants 

I± = u ± 2^/p - aT (5) 



so that Eqs.Q read as 



where 



dl± dl± 
dZ> C± dT 



c± = -I ± + + aT 



(6) 



(7) 
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If we set u 



aT + v, so that 
U + I- 



P 



16 



use the accelerated reference frame T' 
simplicity, Eqs.(|6]) read as 



dJ± ^ n dJ ± 



dZ 



dT 







where 



J± — v ± 2^ 



C 4 



3 1 
4 J± + 4 J t 



(8) 



T — aZ' 2 /2, Z' = Z', and drop primes for 



(9) 



(10) 



For a temporally localized input optical waveform such as a chirp-free square pulse, 
i.e., with p(T, Z = 0) = p f° r 1^1 < an d p(T, Z = 0) = otherwise, Eq.(|4]) may be 
analytically solved up to the point Z = Tq/ y/po i n terms of the well-known Ritter dam- 
break solution [T2],|20j. On the other hand, henceforth we shall restrict our attention to 
the propagation of special pre-chirped optical pulses, as discussed in the next section. 



3. Nonlinear Riemann waves 

A particular solution of Eqs.([9]), known as the Riemann wave, can be obtained if we set 
J_ = 0, so that v = 2y/p. Therefore Eqs.([9l reduce to 



or 

— h r — 

dZ dT 







where T = 3J+/4 = 2v/2 = The solution of Eq.(ll) reads as T(Z,T) = 

T (T — TZT), where r (T) = YiZ = 0, T). In practice it is easier to consider a given 
intial power profile for the the optical pulse. From the first of Eqs.Q, one also obtains 
for the pulse power the following equation 



dz +6 VP dT 



o 

The solution of Eq.(|12j) corresponding to the initial power profile p(Z 
PqP(T /Tq) can be written in implicit form as [21] 



(12) 



0,T) 



p(Z,T)=p P 



T - 3X v /p(Z,T)) /T 



(13) 



or, in more compact notation, as 

where we defined p = pj p , r = T/T , a 



2,/2p /T , and C 



(14) 

Z/^W Q . In the 



following subsection, we will apply Eq.(14) to exactly solve some specific cases of interest 



in nonlinear optics, namely a pulse with an initial power profile of either parabolic, 
Gaussian or hyperbolic secant shape. 
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Figure 1. Contour plot of pulse power vs. fiber length by solving Eq.Q with an 
initial parabolic power profile. 

Consider first the interesting case of a pulse with initial (i.e., at £ = 0) parabolic 
power profile of finite duration 2T (or parabolic cap), that is 

1 — t 2 for \t\ < 1 

/,(-.<, =[)) = { , (15) 

0/or |r| > 1 



One obtains from Eq.(|14|) the solution at any £ 

/A 



r b 



A 



for r > 0, A > 



for t > 0, A < 
£ - VA 



.4 



/or r < 0, A > 



for t < 0, A < 



where 



' A = (l + a 2 ^ 



5 



1 - r 2 ) (l + a 2 C 2 ) - 2a 2 C V 



A = B — AC 



(16) 



(17) 
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Figure 2. Pulse power temporal profiles at the output of different DCF lengths: 
blue solid thick curves indicate numerical solution of the NLSE; green solid thin curves 
show the exact solution of the NSWEs; the red dashed curve shows the input parabolic 
power profile. 

The corresponding temporal chirp profile is always provided by the relationship 



v = lyfp. It is interesting to compare the analytical solution (16) of the NSWEs 
Q with the numerical solution of the original NLSE Q. Let us consider the 
specific case of a dispersion-compensating fiber (DCF) with normal GVD coefficient 
D = — 100ps/(nm- km) (or = 127 ps 2 / km at 1550 ram), and the nonlinear coefficient 
7 = 3.2W- 1 kmr 1 in Eq.Q. 

In Figjl] we illustrate the propagation of a nonlinear Riemann wave with initial 
parabolic power profile by showing the contour plot of the pulse power as a function of 
the distance along the DCF, as it is computed by solving the orig inal NLSE Q. Here 
the input peak pulse power is equal to P = 2 W, and the parabolic pulse full-width 
at half maximum is Tj w h m = 200 ps. Figjl] shows that the pre-chirp distribution leads 
to a slowing down of the pulse. Moreover, the pulse power profile develops a strongly 
asymmetric shape: its trailing edge acquires a vertical slope at about z = 1 km. For 
longer distances both trailing and leading edges of the pulse exhibit a nearly linear 
decrease with time, albeit with different slopes. The details of the evolution of the 
pulse power profile with fiber length is clearly illustrated in Fig|2j here we compare the 
numerical solution of the NLSE (jlj (blue solid thick curves) with the exact solution of 
the NSWEs (green solid thin curves). In Figj2]we also show the input parabolic power 
profile (red dashed curves). 

As it can be seen in Figj2j there is an excellent agreement between the numerical 



solution of the NLSE and the exact solution (16) up to the point (i.e., z = 1 km) where 
a wavefront with vertical slope develops at the trailing edge of the pulse. For longer 
distances, the agreement between the analytical NSWE solution and the numerical 
NLSE solution remains for the entire leading edge of the pulse. On the other hand 
the vertical trailing edge is preserved in the analytical solution of the NSWEs, whereas 
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the trailing edge shows a smoothened decay in the solution of the NLSE. The optical 
pulse deformation which is seen in Figsjl]j2] is analogous to the shoaling of a shallow 
water wave when it approaches to the beach. Therefore we may name the point where 
the vertical pulse edge develops as the shoaling point. It is remarkable that no wave 
breaking occurs in spite of the occurrence of a vertical slope in the pulse trailing edge: 
the NLSE solutions exhibit a self-regularisation behavior. 
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Figure 3. Pulse chirp temporal profiles at the output of different DCF lengths: blue 
solid thick curves indicate numerical solution of the NLSE; green solid thin curves show 
the exact solution of the NSWEs; red dashed curves show the input parabolic power 
profile. 



In Fig|3]we display the evolution of the chirp profiles with distance, corresponding 
to the power profiles of FigfJ} Here the red dashed curves indicate the initial chirp 
profile of the parabolic Riemann pulse. As can be seen, the analytical solutions for 
the chirp, as they are obtained either from the NSWE or the NLSE, are in excellent 
agreement with each other up the shoaling point. For longer distances, the agreement 
only remains for the (relatively longer) leading edge of the pulse. 



3.2. Gaussian pulse 

In the case of a nonlinear Riemann wave with an input Gaussian power profile, we may 



set p{r,( = 0) = exp (— 2r 2 ), so that from Eq.(14) we obtain the implicit equation for 
the pulse power 

p(r, C) = exp -2 (r 2 - 2a(r^ + a 2 C, 2 p) (18) 

which can be easily solved by an iterative procedure. In Figj4] we display the 
corresponding contour plot of the pulse power vs. distance along the DCF, as it is 
obtained by solving the NLSE (jl]). The input Gaussian peak pulse power is always 
equal to P = 2 W, and the pulse Tf whm = 200 ps. Fig|4] shows that with a Gaussian 
pulse the shoaling point occurs at about z = 2.5 km. 
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Figure 4. Contour plot of pulse power vs. fiber length by solving Eq.([T]) with initial 
Gaussian power profile. 




The comparison among the exact NSWE and numerical NLSE solutions for the 
pulse power profiles shown in Fig|5] reveals that their agreement is very good up to the 
shoaling point, that is z = 3 km. For longer distances, the NLSE solution develops a 
tail in the trailing edge which is not captured by the exact solution (18). 

Moreover the comparison of the chirp profiles in Figj6] shows that past the shoaling 
point the NLSE pulse acquires a temporally oscillating negative chirp in its leading 
edge, and a uniform positive chirp in its trailing edge. Such highly chirped regions in 
the pulse tails are not fully captured by the NSWE solution (18). 



Optical fiber tsunamis 



10 




200 

6 100 z = 3 km 



-500 -250 

200 1 — i r- 

100^ z=4km 


-iooL-i ' 

Time [ps] 

Figure 6. Same as in Fig|3j with initial Gaussian pulse power profile. 
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Figure 7. Contour plot of pulse power vs. fiber length by solving Eq.fJTJ) with initial 
hyperbolic secant power profile. 



3.3. Hyperbolic secant pulse 

In the case of an input hyperbolic secant pulse, we may set p(r,( = 0) = sech 2 {j) so 



that, from Eq.(14) we obtain 



p(r,C) = sech 2 (r - aC^/p) (19) 

In Fig|7]we illustrate the contour plot of the pulse power vs. distance along the DCF 
from the solution of the NLSE pi) for an input Riemann wave with initial hyperbolic 
secant power profile. The input peak pulse power is kept equal to P = 2 W, and 
Tfwhm = 200 ps. FigJT] shows that with a hyperbolic secant pulse the shoaling point 
is shifted further up to about z = 3.5 km: otherwise the power evolution remains 
qualitatively very similar to the case of a Gaussian pulse which is shown in Fig|4j 
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Figure 8. Same as in Fig[2j with initial hyperbolic secant pulse power profile. 



Figj8] illustrates the details of the pulse power with an initial hyperbolic secant 
profile for the exact NSWE and numerical NLSE solutions. Again, there is an excellent 
agreement up to the shoaling point, or z = 4 km. 
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Figure 9. Same as in Fig|3j with initial hyperbolic secant pulse power profile. 

On the other hand, the comparison of the chirp profiles which is illustrated in 
Figj9] shows that past the shoaling point the NLSE pulse acquires a strong, temporally 
oscillating positive chirp in its trailing edge. Such highly chirped region corresponds to 



the tail which is not reproduced by the NSWE solution (19). 
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4. Rogue waves 

In the previous section [3] we have shown that properly pre-chirped optical pulses with 
different power profiles experience the wave-breaking free formation of a vertical edge, 
which is analogous to the shoaling of shallow water waves. On the other hand, in 
spite of the strong temporal deformation, the value of the peak power remained nearly 
unchanged and close to the input value in the course of propagation. In this section 
we will demonstrate by numerical simulations that extreme wave formation, that is the 
emergence of temporally compressed, high intensity and transient pulses, is possible in 
optical fibers with the inclusion of frequency shifting and TOD. This is equivalent to 
the case of a dispersion varying optical fiber, where the local value of the dispersion 
changes owing to the variation of the fiber diameter [18j-[19j. Let us consider next the 
impact of fiber TOD on the propagation of nonlinear Riemann waves in optical fibers: 
the corresponding modified NSWE reads as [22] 



dp 
d~Z> 

du 



+ (l + f3u)u^-+ (l + 2(3u)p 



where (3 = /3/(2/3 2 



(l + pu) 



dT 

du dp 

v 1 — 

dT dT 



du 
'dT 







a. 



(20) 
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Figure 10. Same as in Figj6j with negative third-order dispersion 



In the following, we restrict our attention to investigating the possibility of the 
emergence of spatio-temporally localized flashes or rogue-wave type of solutions under 
the action of nonlinearity, dispersion, frequency-shifting and TOD by numerically 
solving the original NLSE ([!]). For simplicity, we will only consider the case of an 
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initial Riemann pulse with a Gaussian power profile, for the same power and temporal 
duration values as in the previous section. 

Fig. 10 shows the contour plot of the pulse power evolution along the DCF whenever 
the frequency up-shift rate is equal atf = 256 MHz/m and the TOD equal to 
/3 3 = —2ps 3 /m. Note that the zero-dispersion point is up-shifted by 31.4 GHz with 
respect to the carrier frequency of the input pulse. The direction of frequency shifting is 
such that the fiber GVD is reduced by TOD as the pulse propagates in the fiber, so that 
the zero- dispersion point is crossed and eventually the pulse moves into the anomalous 
dispersion regime. 

In FigjTU] and in the rest of this section, we will display the pulse evolution in the 
original temporal frame and not in the accelerated frame as it was done in the previous 
section. As it can be seen by comparing Fig. 10 with Fig|4| the relatively large TOD value 
has the effect of anticipating the shoaling point (or vertical trailing edge) to z — 0.3 km 
(from z = 2.5km in the absence of TOD). Moreover, a significant temporal compression 
and three-fold peak power amplification (from 2 W to 6 W) occurs in correspondence of 
the shoaling point. 



7 1 1 1 1 1 1 i 
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Figure 11. Details of pulse evolution as in Fig 11 numbers next to each curve indicate 
the propagation distance in meters. 

A selection of pulse power profiles at a given set of distances (indicated in meters 



next to each curve) along the fiber, corresponding to the power contour plot of Fig 10 
is shown in FigJTTJ Here we can see that the formation of a vertical trailing edge at 
z = 300 m is accompanied by a nearly three-fold amplification of the peak power value. 
Past the shoaling point, the peak power decreases, and a temporal oscillation appears 
in the leading edge while the trailing edge remains smooth. 

Quite suprisingly, numerical simulations reveal that whenever the sign of the TOD 
is changed from negative to positive (i.e., we set (3 3 = 2ps 3 /m), that is the frequency- 
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Figure 12. Same as in Fig 10 with positive third-order dispersion /3 3 = 2ps /m. 



shifting pulse sees a progressively larger normal GVD as it propagates along the fiber, 
the temporal compression and peak power enhancement of the input Gaussian Riemann 
pulse is greatly enhanced at the shoaling point. 
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Figure 13. Surface plot of pulse power vs. distance for the same case as in Fig 12 



In Figs. 12 13 we show the contour and surface plots of the pulse power evolution 



along the DCF for the same frequency up-shift rate as in Figs. 10 but the sign of TOD 
is changed from negative to positive. As it can be seen, in this case the vertical edge 
or shoaling occurs in the leading edge of the pulse at about z = 400 m. Moreover, the 
temporal compression and the pulse peak power is substantially increased with respect 



Optical fiber tsunamis 



15 



to the case of negative dispersion which was illustrated in Figs, 10 11 In particular, a 
more than ten-fold increase in the peak power is observed at the point of maximum 
temporal compression z = 400 m. 



20- 




Figure 14. Details of pulse evolution as in FigfT2| numbers next to each curve indicate 
the propagation distance in meters. 

A set of pulse power profiles for progressively increasing distances (indicated in 



meters next to each curve) in the fiber, extracted from the plots of Fig, 12 13, is shown 



in Fig, 14 As it can be seen, the compressed pulse reaches a peak power as high as 
24 W at z = 400 m. Past the point of maximum compression, the pulse break-up into 
multiple peaks until a pulse train is obtained at z — 600 m. 

In the previous examples we have used a relatively large value of the TOD 
coefficient. However the dynamics of rogue wave formation is not critically dependent 
on the specific TOD value: in Figjl5] we show that a similar pulse compression and 
extreme peak formation also occurs whenever the TOD is reduced by five times to 



/?3 = 0.4 /m. The main difference between Fig 14 and Fig, 15 is that the point of 



maximum compression is moved to further down the fiber, at z = 1 km from z = 400 m. 



Fig|T6] displays the spectral intensities associated with the input pre-chirped pulse 
(red dashed curve) and the frequency up-shifted rogue pulse at z — 1 km which are 
reported in FigfT5j As can be seen, in spite of the relatively large normal GVD, the 
spectrum of the rogue pulse exhibits broad triangular tails that are associated with the 
presence of a sharp hyperbolic secant soliton-like peak. 

It is interesting to point out that, in contrast with the case of either chirp-free 
input pulses, which are subject to TOD-induced wave-breaking or shock formation in 
dispersion tapered fibers [H]-[T9], in the case of Riemann pulses rogue wave formation 
is not accompanied by any wave-breaking phenomena. Note that whenever a step-wise 
pre-chirp of a CW is used, temporal shocks still arise even in the absence of TOD 
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Figure 15. Same as in Fig 14 with the positive third-order dispersion /?3 = OAps /m. 
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Figure 16. Spectra of input pulse and of maximum compressed pulse as in Fig 15 



ng-ng. 

In order to highlight the nonlinear and dispersive mechanisms of pulse formation 
dynamics, in Eq.Q and elsewhere in this paper we neglected the presence of linear 
optical fiber loss. Indeed the effect of fiber loss is minimal and it does not qualitatively 
change our conclusions for the relatively short propagation lengths (up to a few km) 
which have been considered in our numerical examples. In general, effective transparent 
propagation can be achieved in optical fibers thanks to distributed Raman amplification, 
or by means of periodic lumped amplification in optical communication systems. 
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5. Conclusions 

In this work we have described optical pulse shoaling in the normal dispersion regime 
of optical fibers. We obtained exact solutions of the optical NSWE, and demonstrated 
their relatively good match with the numerical solutions of the NLSE. We have also 
revealed that TOD may lead to the occurrence of extreme waves or optical tsunamis 
in dispersion tapered fibers, in full analogy with the dramatic run-up and wave height 
amplification of ocean tsunamis. 

The present results may also have applications in a context different from nonlinear 
optics, such as hydrodynamics or BEC. Moreover, we envisage that optical Riemann 
waves and tsunamis may also be observed in the spatial domain, that is when diffraction 
replaces dispersion. In this case the presence of two-dimensional degrees of freedom may 
facilitate the generation and control of the necessary initial phase profiles [23] . 
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